Nuclear Stopping in Au+Au Collisions at A /s/v r /v = 200 GeV 
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Transverse momentum spectra and rapidity densities, dN/dy, of protons, anti-protons, and net- 
protons (p — p) from central (0-5%) Au+Au collisions at ^snn = 200 GeV were measured with 
the BRAHMS experiment within the rapidity range < y < 3. The proton and anti-proton dN/dy 
decrease from mid-rapidity to y = 3. The net-proton yield is roughly constant for y < 1 at 
dN/dy ~ 7, and increases to dN/dy ~ 12 at y ~ 3. The data show that collisions at this energy 
exhibit a high degree of transparency and that the linear scaling of rapidity loss with rapidity 
observed at lower energies is broken. The energy loss per participant nucleon is estimated to be 
73 ± 6 GeV. 

PACS numbers: 25.75 Dw. 



The energy loss of colliding nuclei is a fundamental 
quantity determining the energy available for particle 
production (excitation) in heavy ion collisions. This 
deposited energy is essential for the possible formation 
of a deconfined quark-gluon phase of matter (QGP). 
Because baryon number is conserved, and rapidity dis- 
tributions are only slightly affected by rescattering in 
late stages of the collision, the measured net-baryon 
(B — B) distribution retains information about the 
energy loss and allows the degree of nuclear stopping 
to be determined. Such measurements can also dis- 
tinguish between different proposed phenomenological 
mechanisms of initial coherent multiple interactions 
and baryon transport Ql 0- 01 ■ 

The average rapidity loss, (Sy) = y p — (y) 1231ms 
used to quantify stopping in heavy ion collisions |j,|5j]. 
Here, y p is rapidity of the incoming projectile and (y) 
is the mean net-baryon rapidity after the collision : 
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where N part is the number of participating nucleons 
in the collision. The two extremes correspond to full 
stopping, where initial baryons lose all kinetic en- 
ergy ((Sy) = y p ) and full transparency, where they 
lose no kinetic energy ((Sy) = 0). For fixed colli- 
sion geometry (system size and centrality) at lower 
energy (SIS, AGS, and SPS) it was observed that (Sy) 
is proportional to the projectile rapidity. For cen- 
tral collisions between heavy nuclei (Pb, Au), (Sy) ~ 
0.58-^0,10. 

Bjorken assumed that sufficiently high energy colli- 
sions are "transparent" , thus the mid-rapidity region 
is approximately net-baryon free . The energy den- 
sity early in the collision, e, can then be related in a 
simple way to the final particle production. At RHIC 
it has been estimated that e ~ 5 GeV/fm 3 , well above 
the lattice QCD prediction (e cr a ~ lGeV/fm 3 [9j) for 
the hadron gas to QGP phase transition. 

In this letter, results on proton and anti-proton 
production, and baryon stopping in Au + Au colli- 
sions at y?SNN = 200 GeV are presented. The data 
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were collected with the BRAHMS detector at RHIC. 
The BRAHMS experiment consists of two indepen- 
dent spectrometer arms, the Mid-Rapidity Spectrom- 
eter (MRS) and the Forward Spectrometer (FS), de- 
scribed in detail in [l^}. The spectrometers consist 
of dipole magnets, Time Projection Chambers (TPC) 
and Drift Chambers (DC) for tracking charged par- 
ticles, and detectors for particle identification (PID). 
The MRS can be rotated 30° < 6 < 95° and the FS 
2.3° < 8 < 30°, where 9 is the polar angle with respect 
to the beam axis. By combining different settings 
of angle and magnetic fields, (anti-)proton transverse 
momentum spectra at different rapidities (0 < y < 3) 
were obtained. 

The interaction point (IP) is determined by the 
Beam-Beam counters (BB), with a precision of <jbb = 
0.7 cm. The IP was required to be within 15 (20) cm 
of the nominal IP for the MRS (FS) analysis to mini- 
mize acceptance corrections. 

Protons and anti-protons are identified using Time 
Of Flight detectors (TOF) in the MRS (TOFW) and 
FS (HI and H2), and by the Ring Imaging Cherenkov 
(RICH) in the FS. The TOF resolution is a TO FW = 
80 ps, am = &H2 — 90 ps. To identify (anti-)protons 
it is required that the derived square of the mass, m 2 , 
is within ±2er of the proton m 2 and at least 2a away 
from the m 2 for kaons. This allows (anti-) protons 
to be identified up to momentum p < 3.0 GeV/c for 
TOFW, p < 4.5 GeV/c for HI, and p < 6.5 GeV/c 
for H2. The refractive index of the RICH, n = 1.0020 
allows protons to be directly identified via a measure- 
ment of their ring radius in the range 15 < p < 25 
GeV/c. By using the RICH to veto pions and kaons 
the PID can be extended down to p = 10 GeV/c. 

The collision centrality is determined using a mul- 
tiplicity array located around the nominal IP [10| . In 
this analysis only the 0-5% most central collisions were 
used. The mean number of participants for this cen- 
trality class was found by Glauber model calculations 
to be N part = 357 ±8 [13 • 

From the identified particles, invariant differential 
yields 55^7 dydp T were constructed for each spectrom- 
eter setting, e differential yields were corrected for 
geometrical acceptance, tracking and PID efficiency, 
absorption and multiple scattering. The corrections 
for absorption and multiple scattering were less than 
20% at the lowest pr and less than 5% at the high- 
est pr in all settings. No corrections were applied for 
secondary protons from e.g. the beam pipe, since the 
contribution was found to be negligible in GEANT 
based Monte Carlo (MC) simulations 0] using HI- 
JING as input, when the tracks were required 
to point back to the IP. The acceptance correction is 
purely geometrical and calculated using a MC sim- 
ulation of the BRAHMS detector. The tracking ef- 
ficiency was estimated for the TPCs, by embedding 



simulated tracks in real events (method I) , and in the 
FS where there are 5 tracking detectors (2 TPCs and 
3 DCs) by comparing the number of identified track 
segments in the chamber to the number of reference 
tracks determined by other detectors disregarding the 
chamber under consideration (method II). In the front 
part of the FS (FFS) where both methods are appli- 
cable, there was a 10% discrepancy between the two 
methods which has been included in the systematic 
errors quoted below. For the spectrometers the total 
tracking efficiency depends on the spectrometer angle 
and is 90-95% in the MRS, 80-90% in the FFS (set- 
tings at y ~ 2) and 60-70% in the full FS {y ~ 2, 
y ~ 3). The efficiencies of the TOF detectors were 
found to be 93-98%. For the part of the spectrum 
at y ~ 3, where the RICH was used to veto pions 
and kaons a correction for contamination of pions and 
kaons with no identified ring radius was applied. The 
RICH inefficiency is 5%. 
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FIG. 1: Proton, anti-proton, and net-proton transverse 
momentum spectra at selected rapidities, y ~ 0,1,2,3. 
The solid and dashed lines represents Gaussian fits to the 
data (extended outside the fit-range for clarity). The anti- 
proton spectra have been divided by 10 as indicated. The 
error on the data points are statistical only and no weak 
decay correction has been applied. 

Figure ^ shows transverse momentum spectra for 
four of the nine measured rapidities. The net-proton 
spectra were constructed by subtracting the anti- 
proton spectrum from the proton spectrum. At each 
rapidity, the proton, anti-proton, and net-proton 
spectra have similar shape, indicating that produced 
and transported protons have similar spectral proper- 
ties. To obtain rapidity densities, dN/dy, for p,p and 
net-p their spectra are fitted, and the fit was used to 
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extrapolate to the full pr range. Different functional 
forms were tested: my-exponential, Boltzmann and 
Gaussian. The function found to best describe the 
data was the Gaussian in p? [f(p T ) oc e -P5V(20] a nd 
this function has been used for all fits. This functional 
form was also used in (lif . 

The mean transverse momentum (pr) of the spectra 
calculated from the fit is found to be within 0.1 GeV/c 
at each rapidity for the three functional forms. For 
protons which have the best counting statistics, (pr) 
decreases from (pr) = 1.01 ± O.Ol(stat) GeV/c at y = 
to {p T ) = 0.84 ± 0.01 (stat) GeV/c at y ~ 3. 
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units of rapidity, from dN/dy(y=0) = 6.4±0.4(stat)± 
l.O(syst) to dN/dy(y=3) = 12.4±0.3(stat)±3.2(syst). 

A Gaussian fit to the anti-proton dN/dy distribu- 
tion gives the total extrapolated anti-proton yield : 
84±6 (92% in -3 < y < 3). For protons the yield from 
a Gaussian fit to dN/dy in the range, — 3 < y < 3, is 
138 ±7. 
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FIG. 3: The net-proton rapidity distribution at AGS 0, 
DHH (Au+Au at v^7= 5 GeV), SPS £3 (Pb+Pb at 
y / ijvjv = 17 GeV) and this measurement (-^/sjvjv = 200 
GeV). The data are all from the top 5% most central 
collisions and the errors are both statistical and system- 
atic (the light grey band shows the 10% overall normal- 
ization uncertainty on the E802 points, but not the 15% 
for E917). The data have been symmetrized. For RHIC 
data black points are measured and grey points are sym- 
metrized, while the opposite is true for AGS and SPS data 
(for clarity) . At AGS weak decay corrections are negligible 
and at SPS they have been applied. 



FIG. 2: Proton, anti-proton, and net-proton rapidity den- 
sities dN/dy as a function of rapidity at y / ijvjv = 200 
GeV. The horizontal bars shows the rapidity intervals for 
the projections. The errors shown with vertical lines are 
statistical only while the caps includes both statistical and 
systematic. No weak decay correction has been applied. 

The differential yield within the measured pr range 
varies from 85% of the total dN/dy near mid-rapidity 
to 45% at y ~ 3. The systematic errors on dN/dy 
were estimated from the difference in dN/dy values 
obtained using different spectrometer settings cover- 
ing the same (y,Pr) regions, the discrepancy between 
the two different efficiency methods, and by estimat- 
ing the effects of the pr extrapolation. The system- 
atic errors were found to be 10-15% for mid-rapidity 
(y < 1) and 20-30% for forward rapidities. 

Figure |2 shows the resulting rapidity densities 
dN/dy as a function of rapidity. The most prominent 
feature of the data is that while the proton and anti- 
proton dN/dy decrease at rapidities away from mid- 
rapidity the net-proton dN/dy increases over all three 



Figurc|3shows net-proton dN / dy measured at AGS 
and SPS compared to these results. The distributions 
show a strong energy dependence, the net-protons 
peak at mid-rapidity at AGS, while at SPS a dip is 
observed in the middle of the distribution. At RHIC a 
broad minimum has developed spanning several units 
of rapidity, indicating that at RHIC energies collisions 
are quite transparent. 

To calculate the rapidity loss, dN/dy must be 
known from mid-rapidity to projectile rapidity, y p — 
5.36. BRAHMS measures to y ~ 3, so the shape of the 
rapidity distribution must be extrapolated to calcu- 
late (5y). The baryon number of participating nucle- 
ons (Np ar t) is conserved, while the net-proton number 
is not necessarily conserved. To obtain net-baryons, 
the number of net-neutrons and net-hyperons have 
to be estimated and the contribution from weak de- 
cays included in the measured net-protons has to be 
deduced. Using MC simulations we find these contri- 
butions to be Ci = 0.53 ±0.05 protons for each A, and 
C2 = 0.49 ± 0.05 protons for each E + decay. There is 
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a weak rapidity dependence which is included in the 
systematic error. The number of net-baryons is, then, 



n B 



nA + «E+ ± 



b p, meets 



■ Cin A + c 2 n s + 



(2) 



where n x — N x — N s are the primary number of 
net-neutrons (n n ), net-protons(n p ), net-lambdas (n\), 
and net-sigmas(ns- , n-£+), respectively, and n p>meas 
is the measured net-proton yield. In addition to pri- 
maries, net-lambdas, n\ includes contributions from 
other hyperons that decay to protons through As e.g. 
£°,S°,andS-. 

The ratio n n /n p = 1.00 ± 0.05 was found from 
HIJING 13] and AMPT 18] in the rapidity inter- 
val \y\ < 3.5. The equilibration of protons and neu- 
trons close to mid-rapidity has been experimentally 



observed at AGS energies 
At mid-rapidity the ratio 
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N A /N p = 0.89 ± 



0.07(stat) ± 0.21(syst) was found to be equal within 
statistical errors to N A /N p = 0.95 ± 0.09(stat) ± 
0.22(syst) at ^/sJTn = 130 GeV Those results 

indicate that N\/N p = N^/Np = n^jn p . We use 
n A /n p = 0.93±0.11(stat)±0.25(syst) at y/JJJ^ = 200 
GeV and find at y — the number of feed-down cor- 
rected protons and anti-protons to be dN/dy — 17. 5± 
1.2(stat) ± 3.0(syst) and dN/dy = 13.2 ± 0.9(stat) ± 
2.3(syst) respectively, in agreement with the measure- 
ment at y/s^N = 200 GeV by PHENIX : dN/dy(p) = 
18.4±2.6(syst),d/V/<2y(p) = 13.5 ± 1.8(syst) 21]. As- 
suming that n\/n p is constant over the full rapidity 
interval the feeddown correction for A alone can be 
done cLS Tip — 0.67 ±0.12 

■ n p,meas- At forward rapid- 
ity As have not been measured, and the conservative 
estimate n\/n p — 1.0 ± 0.5 is used. 

The yield of S + and £~ have not been mea- 
sured at RHIC. In thermal models N%+ ~ and 
iV s - ~ 0.1 -N p H3 at RHIC (^/sWn = 130 GeV), and 
n^+/n p — n-£-/n p — 0.10 ±0.05 is used here. 

The final correction is tib = 2.03 ± 0.08 • n p ^ meas . 
The systematic error on the final correction depends 
almost entirely on the error on n n /n p and n s -/n p 
because the partial derivative of Eq. |21 with respect 
to n\/n p and n s + jn v is very small for the numerical 
values used. 

Figure 0| (insert) shows the net-baryon dN/dy ob- 
tained from the measured net-proton dN/dy. The 
net-baryon density in the rapidity region where 
BRAHMS has no acceptance can be constrained by 
the total integral, equal to N part , and since diffrac- 
tive scatterings can be ignored in central collisions 
where each nuclcon on the average experiences sev- 
eral scatterings, dN/dy (y p ) ~ 0. To illustrate the ex- 
trapolation from the measured region, < y < 3, to 
beam rapidity, y p , two different functions have been 
used: A six order symmetric polynomial (pol6) f(y) 
which is the simplest polynomial that describes the 
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FIG. 4: The inserted plot shows the extrapolated net- 
baryon distribution (data points) with fits (represented 
by the curves) to the data, see text for details. The full 
figure shows the rapidity loss, obtained using Eq. as a 
function of projectile rapidity (in the CM). The hatched 
area indicates the unphysical region and the dashed line 
shows the phenomenological scaling (Sy) = 0.58 ■ y p . The 
data from lower energy are from 0,L|- 



data points, has the correct integral, and f(y p ) = 0, 
and a Bjorken inspired symmetric sum of 2 Gaus- 
sians in momentum space (momgaus), converted to 
rapidity space using p = m sinh(y) (j>t = 0) where 
m is the proton mass. The momgaus gives a sim- 
ple description of the data as the product of the two 
fragmenting nuclei with no mid-rapidity source. The 
data are consistent with Bjorken's picture of the col- 
lision where the fragmentation regions are far from 
mid-rapidity (y ~ 3-4), but shows that there is still 
a significant number of baryons transported to mid- 
rapidity. From the momgaus fit a mean momentum 
of 31 GeV/c with a spread of 13 GeV/c is found. The 
rapidity loss is (Sy) — 2.04 ± 0.10 using the pol6 fit 
and (Sy) = 2.06 ± 0.16 for the momgaus fit, and the 
last value has been used as estimate in Fig. 0] From 
the pol6 fit there are 60 ±6 net-baryons in the rapidity 
interval, < y < 3, with a (Sy) = 3.60±0.05. If the re- 
maining 118.5 net-baryons are placed at y=3.5, maxi- 
mum (Sy) = 2.45, and y — 5.0, minimum (Sy) — 1.45, 
a conservative estimate for the possible range of values 
is obtained, also shown in Fig.^J The rapidity loss at 
RHIC is clearly less than the phenomenological linear 
rapidity scaling would predict; this scaling is broken 
at RHIC. 

Using the functional forms the total energy, E, per 
net-baryon after the collision can be derived: 



E = 
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The (tot) is known for protons from the spectra in 
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the covered rapidity interval. Linear and Gaussian 
extrapolations of (mj) to projectile rapidity changes 
E by less than 5%. The energy is E = 30 ± 2 GeV 
for pol6 and E = 26 ± 5 GeV for the momgaus fits. 
Taking E = 27±6, AE = 73±6 GeV of the initial 100 
GeV per participant is available for excitations. Using 
the same minimum and maximum estimates above : 
47 < AE < 85 GeV. 

In conclusion BRAHMS has measured proton, anti- 
proton, and net-proton yields from mid-rapidity (y = 
0) to forward rapidity (y ~ 3). The net-proton dis- 
tribution shows that collisions at RHIC energies are 
quite transparent compared to lower energies. By ex- 
trapolation to the full net-baryon distribution, we find 
that the rapidity loss scaling observed at lower energy 
is broken and the rapidity loss seems to saturate be- 
tween SPS and RHIC energies. 
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